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1 NZOIROWBI^D INTERFERENCE EXCISION DEVICE 

2 BACKGROUND OF THE mVENTION 

3 The present invention is generally in the field of wireless communication. More 

4 specifically, the invention is in the field of wireless communication receivers using interference 

5 removal. 

6 Wireless communication systems include transmitters capable of transmitting radio 

7 frequency (rf) signals to receivers. Specifically, transmitters transmit a desired rf signal and 

8 receivers receive a received rf signal comprising the desired rf signal, noise and interference. 

9 Noise is defined as non-local (with respect to the receiver) unwanted rf signals, which can be 

10 caused by, for example, biological sonar, rf jammers and solar flares. Interference is defined as 

11 local (with respect to the receiver) unwanted rf signals, which can be caused by, for example, 

12 electromagnetic radiation from electrical equipment and power lines in close proximity to 

13 receivers. Wireless communication systems attempt to improve reception quality by eliminating 

14 or reducing interference. 

15 Hardware notch filtering is a typical method for reducing narrowband interference in 

1 6 wireless communication system receivers. Hardware notch filtering uses an external reference 

17 signal to create notch filters notched at the fundamental and harmonic frequencies of the extemal 

18 reference signal. The hardware notch filters operate in the time domain to remove or reduce 

19 narrowband interference from received rf signals. Disadvantageously, hardware notch filtering 
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1 requires external references, cannot remove strong narrowband interference and can require 

2 multiple notch filters. 

3 Typical excision devices attempt to reduce narrowband interference in wireless 

4 communication system receivers by operating in the frequency domain, which does not require 

5 external references. Typical excision devices include a transformer, a power estimator, an 

6 excisor and an inverse transformer. The transformer performs "windowing" (e.g., Blackman 

7 windowing) and the well-known Fast-Fourier Transform (FFT) on a received signal to transform 

8 the received signal from the time domain to the fi-equency domain. Windowing techniques are 

9 described in more detail in an article pubUshed by the Proceedings of the IEEE in 1978, entitled 

1 0 "On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform", authored 

11 by Frederic J. Harris, which is hereby incorporated by reference in its entirety for its teachings on 

12 windowing techniques. The FFT is described in more detail in a book published by Macmillan 

13 in 1988, entitled "Introduction to Digital Signal Processing", authored by J. G. Proakis et al, 

14 which is hereby incorporated by reference in its entirety for its teachings on the FFT. The power 

15 estimator receives fi-equency domain data from the transformer and determines a power vector. 

1 6 The excisor receives fi-equency domain data from the transformer and the power vector from the 

17 power estimator. The excisor removes narrowband interference by determining whether the 

18 power vector exceeds an excision threshold. The inverse transformer receives frequency domain 

19 data from the excisor and inverse transforms the data into time domain data. Disadvantageously, 
2 0 typical excision devices cannot adapt to noisy changes in environment. For example, high 

2 1 impulsive noise (i.e., short duration, high-ampUtude noise) and high interference precludes 

22 typical excision devices from removing narrowband interference in Extremely Low Frequency 
2 3 (ELF) communications. 

24 Therefore, a need exists for narrowband interference excision devices that do not require 

25 an external reference for narrowband interference. In addition, a need exists for a method and 

2 6 apparatus for narrowband interference excision devices that operates in high impulsive noise and 

2 7 high interference environments. 

28 
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1 SUMMARY OF THE INVENTION 

2 The present invention is directed to a method and apparatus for narrowband interference 

3 excision. The invention overcomes the need in the art for narrowband interference excision 

4 devices that do not require an external reference for narrowband interference and that operate in 

5 high impulsive noise and high interference environments. 

6 The present inventive narrowband interference excision device method and apparatus 

7 uses a transient eUminator and a moving average calculator to remove high-amplitude 

8 interference. The present invention is capable of removing high-amplitude interference in high 

9 impulsive noise and high-interference environments without external references. The method 

10 and apparatus is particularly useful in ELF communication systems 

11 According to one embodiment, the present invention is a method for narrowband 

12 interference excision. The method comprises a step of transforming time domain data into a 

13 frequency domain vector having a plurality of frequency bins. Next, the method estimates a 

14 plurality of power values corresponding to frequency bins of the frequency domain vector to 

15 obtain a current power vector. Then, the method calculates an average power vector from the 

1 6 current power vector and at least one selected previous power vector. Further , the method 

17 excises selected frequency bins of the fi-equency domain vector to produce an excised frequency 

1 8 domain vector. Moreover, the method determines whether to include the average power vector 

19 in fiiture average power vector determinations. 

2 0 According to another embodiment, the present invention is a narrowband interference 

2 1 excision device, where the narrowband interference excision device includes a transformer, 

2 2 power estimator, moving average calculator, excisor and transient eliminator. The transformer ii 

2 3 capable of transforming time domain data to a frequency domain vector comprising a plurality o 

2 4 frequency bins. The power estimator is operatively coupled to the transformer and capable of 

2 5 receivmg the frequency domain vector and adapted to determine a PXW power vector 

2 6 comprising a pluraUty of power values that correspond to the plurality of frequency bins. The 

2 7 moving average calculator is operatively coupled to the power estimator and is capable of 

2 8 receiving the PXW power vector and a discard flag. The moving average calculator is adapted 1 

2 9 determine an average power vector. The excisor is operatively coupled to the transformer and h 

3 0 capable of receiving the frequency domain vector and the average power vector. The excisor is 
3 1 adapted to excise selected frequency bins of the plurality of frequency bins and determine an 

3 
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excised power vector and an excision flag vector. The transient eliminator is operatively coupled 

2 to the excisor and the moving average calculator and is capable of receiving the excision flag 

3 vector. The transient eliminator is adapted to determine whether to include the average power 
vector from fixture average power vector determinations and determine the discard flag. 

The previously summarized features and advantages along with other aspects of the 
present invention will become clearer upon review of the following specification taken together 

7 with the included drawings. 

8 RRIRF DESCPIPTION OF TH K DRAWING 

9 FIG. 1 is a block diagram representing a narrowband interference excision device 

1 0 according to one embodiment of the invention. 

11 FIG. 2 is a block diagram representing a transformer of a narrowband interference 

1 2 excision device according to one embodiment of the invention. 

13 FIG. 3 is a block diagram representing a moving average calculator of a narrowband 

14 interference excision device according to one embodiment of the invention. 

15 FIG. 4 is a block diagram representing a transient eliminator of a narrowband interference 

1 6 excision device according to one embodiment of the invention. 

17 FIG. 5 is a block diagram representing an inverse transformer of a narrowband 

1 8 interference excision device according to one embodiment of the invention. 

1 9 FIG. 6 is a flowchart of an exemplary method that implements an embodiment of the 

20 invention. 

2 1 FIG. 7 is a flowchart of an exemplary method that implements an embodiment of the 

22 invention. 

2 3 FIG. 8 is a flowchart of an exemplary method that implements an embodiment of the 

24 invention. 

2 5 FIG. 9 is a flowchart of an exemplary method that implements an embodiment of the 

26 invention. 

27 FIG. 1 0 is a block diagram representing an exemplary ELF receiver application of a 

28 narrowband interference excision device according to one embodiment of the invention. 
29 



4 



Express Mail No.: EE676641076US 



Navy Case No. 84436 



1 DETAILED DESCRIPTION OF THE INVENTION 

2 The present invention is directed to a method and apparatus for narrowband interference 

3 excision devices. Although the invention is described with respect to specific embodiments, the 

4 principles of the invention, as defined by the claims appended herein, can obviously be applied 

5 beyond the specifically described embodiments of the invention described herein. Moreover, in 

6 the description of the present invention, certain details have been left out in order to not obscure 

7 the inventive aspects of the invention. The details left out are within the knowledge of a person 

8 of ordinary skill in the art. 

9 The drawings in the present application and their accompanying detailed description are 

1 0 directed to merely exemplary embodiments of the invention. To maintain brevity, other 

11 embodiments of the invention that use the principles of the present invention are not specifically 

12 described in the present application and are not specifically illustrated by the present drawings. 

13 The present inventive narrowband interference excision device method and apparatus 

14 uses a transient eliminator and a moving average calculator to remove high-amplitude 

1 5 interference. The present invention is capable of removing high-amplitude interference in high 

1 6 impulsive noise and high-interference environments without external references. The method 

17 and apparatus is particularly useful in ELF communication systems. 

1 8 FIGURE 1 is a block diagram representing a narrowband interference excision device 

19 according to one embodiment of the invention. As shown in FIG. 1, narrowband interference 
2 0 excision device 100 includes transformer 1 10, power estimator 120, moving average calculator 

21 130, excisor 1 40, transient eliminator 150, inverse transformer 1 60 and output processor 1 70. 

22 Transformer 1 10 receives time domain data, assembles time domain data blocks and transforms 
2 3 time domain data blocks into fi-equency domain data. In one embodiment, transformer 1 1 0 uses 
24 a FFT to transform time domain data into fi-equency domain data. Those skilled in the art shall 

2 5 recognize that other methods that transform time domain data into fi-equency domain data can be 

2 6 used with the present invention without departing from the scope and spirit of the present 

2 7 invention. In one embodiment, transformer 1 10 uses an "overlap-and-add process" and a 

2 8 windowing technique to assemble time domain data blocks. Transformer 1 10 inputs frequency 

2 9 domain data blocks to power estimator 1 20 and excisor 140. 

3 0 FIGURE 2 is a block diagram representing a transformer of a narrowband interference 

3 1 excision device according to one embodiment of the invention. As shown in FIG. 2, transformer 
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210 includes current input data block assembler 212, window coefficient generator 214, 



1 

2 multiplier 216 and time-to-frequency transformer 218. Current input data block assembler 212 

3 receives time domain data and assembles the data into data blocks (i.e., vectors) having a length 

4 N. In one embodiment, N is selected as a power of 2 to increase transformation efficiency and 

5 provide adequate frequency resolution and bin signal-to-noise ratio (SNR) for specific 

6 applications. For example, N equals 1 6K (i.e., 1 6* 1 024) in an ELF receiver application. In one 

7 embodiment, assembler 212 uses a "75% overlap-and-add process" to assemble dala blocks from 

8 time domain data. The 75% overlap-and-add process appends 25% of a new data block to 75% 

9 of a previous data block to form data vector "x" having a length N. Assembler 2 1 2 inputs data 

1 0 vector x to multiplier 216. 

11 Window coefficient generator 214 of transformer 210 of FIG. 2 inputs a window 

12 coefficient vector "w" to multipUer 216. In one embodiment, generator 214 uses a Blackman 

13 window coefficient vector according to the following Equation 1 : 
14 

15 (Equation 1) w(k) = 0.42 + (0.50 *cos(2nn/iV))+ (0.08 *cos(4n«/ AO) ; 

16 where <«<-—; 

2 2 

17 w(n) = Blackman window coefficient; and 

18 iV = window size. 
19 

20 In one embodiment, the Blackman window coefficient vector is generated as a periodic function 

2 1 having no symmetry at end points. Generator 2 14 inputs window coefficient vector w to 

22 multiplier 216. 

2 3 Multiplier 2 1 6 receives data vector x from assembler 2 1 2 and window coefficient vector 

24 w from generator 214. Multiplier 2 1 6 multiplies data vector x by window coefficient vector w to 

2 5 obtain windowed vector "xw" and inputs windowed vector xw to time-to-frequency fransformer 

26 218. 

2 7 Transformer 2 1 8 fransforms windowed vector xw from the time domain to the frequency 

2 8 domain to produce frequency domain vector "XW," which comprises a plurality of frequency 

2 9 bins. In one embodiment, fransformer 218 uses a FFT having a block size N according to the 

3 0 following Equation 2: 
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1 

2 (Equation 2) -^(^) = Z ^(«)^^'''""*'''^ ; 

3 where 0<A;<(iV-l) 

4 = frequency domain output; 

5 x(k) = time domain input; and 

6 iV=FFT block size. 

7 

8 Transformer 2 1 8 outputs frequency domain vector XW. 

9 Referring again to FIG. 1 , power estimator 120 receives frequency domain vector XW 

1 0 from transformer 110. Power estimator 1 20 estimates power values of frequency bins of 

11 frequency domain vector XW. In one embodiment, power estimator 120 multiplies frequency 

12 domain vector XW by its complex conjugate to obtain power vector "PXW." Thus, a power 

13 value is computed for each frequency bin. Power vector PXW comprises a plurality of power 

14 values corresponding to frequency bins. Power estimator 120 inputs power vector PXW to 

15 moving average calculator 130. 

16 As shown in FIG. 1 , movmg average calculator 130 receives power vector PXW from 

17 power estimator 120 and discard flag "F" from transient eliminator 150. According to the 

1 8 invention, calculator 1 30 scales power values of frequency bins of power vector PXW using a 

1 9 power scale factor vector PSFV to obtain scaled power vector "PXW," which reduces 

2 0 suppression of high-amplitude frequency bins of the current power vector PXW when calculating 

2 1 average power vector PXW_avg. Calculator 1 30 determines average power vector PXW_avg 

2 2 from scaled power vector PXW, previous power vectors PXW and power scale factor vector 

2 3 PSFV. Calculator 1 30 uses discard flag F to determine whether to discard or retain the current 

2 4 power vector PXW in fiiture calculations of average power vector PXW_avg. Calculator 1 30 
2 5 inputs average power vector PXW_avg to excisor 140. 

2 6 FIGURE 3 is a block diagram representing a moving average calculator of a narrowband 

2 7 interference excision device according to one embodiment of the invention. As shown in FIG. 3, 
2 8 moving average calculator 330 includes power scale factor computer 332, power scaler 334, 

2 9 current block power rescaler 336, previous power block storage 338 and average power vector 

3 0 computer 340. Power scale factor (PSF) computer 332 receives current power vector PXW and 
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1 determines a power scale factor PSF for frequency bins of current power vector PXW. In one 

2 embodiment of PSF computer 332, if the power value of a frequency bin is less than a threshold 

3 "T," the PSF for the frequency bin is set to 1 , else the PSF is set to "M." M is the number of 

4 blocks used for averaging in average power vector computer 340 as described below. In one 

5 embodiment, M ranges between 5 and 10, inclusive. The PSF computer 332 calculates threshold 

6 T for frequency bins. In one embodiment, threshold T is calculated by the following Equation 3: 



7 

8 (Equations) T = MxPXW _avg; 

9 where T - threshold for a frequency bin; 

10 M- number of blocks used for averaging; and 

1 1 PXWjxvg ~ average power for a frequency bin. 
12 



13 Power scale factor vector PSFV comprises power scale factors corresponding to frequency bins. 

14 PSF computer 332 inputs power scale factor vector PSFV and current power vector PXW to 

1 5 power scaler 334. 

1 6 Power scaler 334 receives current power vector PXW and power scale factor vector 

17 PSFV. Power scaler 334 determines a scaled power vector PXW\ In one embodiment, power 

18 scaler 334 receives current power vector PXW from power estimator 120 of FIG. 1. In one 

19 embodiment, power scaler 334 receives current power vector PXW from PSF computer 332. 
2 0 Power scaler 334 receives power scale factor vector PSFV from PSF computer 332. Scaled 

2 1 power vector PXW* is calculated by multiplying power vector PXW by power scale factor vector 

22 PSFV. Power scaler 334 inputs scaled power vector PXW' to average power vector computer 



23 340 and current block power rescaler 336. 

24 Average power vector computer 340 receives scaled power vector PXW' and previous 

2 5 power vectors PXW. Average power vector computer 340 determines an average power vector 

2 6 PXW_avg. In one embodiment, average power vector PXW^avg is calculated by the following 

27 Equation 4: 



28 

-(M-l) 

29 (Equation 4) PXW _avg^ 
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1 where PJffF_avg = average power vector 

2 PXW' = scaled power vector 

3 PXW(n) - n**" power vector 

4 M- number of blocks used for averaging 



5 

6 According to the invention, scaled power vector PXW* helps avoid suppressing high-amplitude 

7 power values in frequency bins of power vector PXW when determining average power vector 

8 PXW_avg. Average power vector computer 340 outputs average power vector PXW_avg. 

9 Current block power rescaler 336 receives scaled power vector PXW' and discard flag F. 

1 0 Current block power rescaler 336 determines power vector PXW. Current block power rescaler 

11 336 can determine power vector PXW by rescaling scaled power vector PXW\ In one 

12 embodiment, power vector PXW is calculated by dividing scaled power vector PXW' by power 

13 scale factor vector PSFV. In one embodiment, current power block rescaler 336 determines 

1 4 power vector PXW by receiving power vector PXW from power estimator 1 20 of FIG. 1 . If 

15 discard flag F indicates retaining current power vector PXW for fiiture averaging by average 

1 6 power vector computer 340, then current power block rescaler 336 inputs current power vector 

17 PXW to previous power block storage 338. 



18 Previous power block storage 338 receives power vectors PXW from current block power 

19 rescaler 336. Storage 338 inputs previous power vectors PXW to average power vector 
2 0 computer 340 for use in determining average power vector PXW_avg. 

21 Referring again to FIG. 1, excisor 140 receives frequency domain vector XW from 



22 transformer 1 10 and average power vector PXW^avg from moving average calculator 130. 

23 Excisor 140 excises selected frequency bins of frequency domain vector XW based on moving 

2 4 averages of frequency bin values of average power vector PXW_avg. Excisor 140 determines an 

2 5 excised frequency domain vector XW' and excision flag vector fExcV, Excisor 140 inputs 

2 6 excised frequency domain vector XW' to inverse transformer 160. Excisor 140 inputs excision 



2 7 flag vector fExcV to transient eliminator 1 50. 

28 In one embodiment of excisor 140, average power vector PXW_avg is used to determine 

2 9 median vector MED, which is used to determine whether to excise selected frequency bins. 

3 0 Median vector MED comprises median values corresponding to frequency bins, where a median 
3 1 value for a specific frequency bin is determined by calculating the median of a set of frequency 
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1 bins that includes the specific frequency bin. For example, the median value for a specific 

2 frequency bin is detennined by calculating the median of a set of frequency bins comprising 

3 seven frequency bins that are centered on the specific frequency bin (i.e., the specific frequency 

4 bin and the three frequency bins immediately lower and three frequency bins immediately higher 

5 than the specific frequency bin). 

6 Excisor 140 compares average power vector PXW_avg to scaled median vector SMED. 

7 Scaled median vector SMED is calculated by multiplying median vector MED by median 

8 threshold MF_T. The value of median threshold MF_T is selected based on the requirements of 

9 an application. Median threshold MF_T should be set high enough to prevent erroneous 

1 0 frequency bin excisions due to random noise. In one embodiment, median threshold MF_T 

11 ranges between approximately 6dB and approximately 12dB, inclusive. 

12 Excisor 140 determines whether an average power value of a frequency bin of average 

1 3 power vector PXW_avg exceeds a scaled median value of a frequency bin of scaled median 

1 4 vector SMED, and if so, excision flag ffixc is set to 1 , otherwise, to 0. A frequency bin of 

1 5 frequency domain vector XW is excised if current excision flag £Exc and previous excision flag 

1 6 fExc are both set to 1 . Excising a frequency bin based on two consecutive excision flags fExc 

17 helps reduce erroneous excisions due to random noise and other transients. Excision flag vector 

18 fExcV comprises current excision flags fExc corresponding to frequency bins. Excision flag 

1 9 vector fExcV is input to transient eliminator 150. Excised frequency domain vector XW 
2 0 comprises unexcised frequency bins of frequency domain vector XW. Excised frequency 

2 1 domain vector XW is input to inverse transformer 160. 

22 As shown in FIG. 1 , transient eliminator 150 receives excision flag vector ffixcV from 
2 3 excisor 140. According to the present invention, transient eliminator 150 determines excision 

2 4 percentage pExc from excision flag vector fExcV and determines whether current power vector 
2 5 PXW should be included in future calculations of average power vector PXW_avg. In one • 
2 6 embodiment, transient eliminator 150 detects transient conditions based on a ratio between 
2 7 current excision percentage pExc and previous excision percentage pExc. When transient 
2 8 conditions are detected, transient eliminator 1 50 sets discard flag F to 1 , which indicates that the 

2 9 current power vector PXW should not be included in fiiture calculations of average power vector 

3 0 PXW_avg. In one embodiment, transient eliminator 1 50 limits the number of power vectors 
3 1 PXW that can be consecutively excised. For example, transient eliminator 1 50 can hmit 
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1 consecutive excision to only three consecutive power vectors PXW. Transient eliminator 150 

2 inputs discard flag F to moving average calculator 130. 

3 FIGURE 4 is a block diagram representing a transient eliminator of a narrowband 

4 interference excision device according to one embodiment of the invention. As shown in FIG. 4, 

5 transient eliminator 450 includes excision percentage calculator 452 and transient excision 

6 device 454. Excision percentage calculator 452 receives excision flag vector fExcV having a 

7 length N. Excision percentage calculator 452 calculates excision percentage pExc by dividing 

8 the summation of excision flags flExc of excision flag vector fExcV by length N. Excision 

9 percentage calculator 452 inputs excision percentage pExc to transient excision device 454. 

10 As shown in FIG. 4, transient excision device 454 receives excision percentage pExc. 

11 Transient excision device 454 determines the ratio between current excision percentage pExc and 

12 previous excision percentage pExc and compares this ratio to a transient condition threshold. If 

13 the ratio is greater than the transient condition threshold, exceed counter cExc is incremented by 

14 1 and transient excision device 454 determines whether to set discard flag F to 1. Discard flag F 

15 is set to 1 if the ratio is greater than the transient condition threshold and exceed counter cExc is 

1 6 less than a predetermined limit such as three. If the ratio is less than or equal to the transient 

17 condition threshold, the exceed counter is set to 0 and current power vector PXW is retained (i.e., 

1 8 used in future average power vector PXW avg calculations). When discard flag F is set to 1 , 

19 moving average calculator 130 excises the current power vector PXW from fiiture average power 
2 0 vector PXW_avg calculations. Transient eliminator 450 inputs discard flag F to moving average 

21 calculator 130 of FIG. 1. 

22 Referring again to FIG. 1, inverse transformer 160 receives excised frequency domain 

23 vector XW' from excisor 140. Inverse transformer 160 transforms excised frequency domain 

24 vector XW' to time domain output vector y'. In one embodiment, inverse transformer 160 uses 

25 an overlap-and-add process to assemble time domain data blocks for output. Inverse transformer 
2 6 160 inputs time domain output vector y to output processor 170. 

2 7 FIGURE 5 is a block diagram representing an inverse transformer of a narrowband 

2 8 interference excision device according to one embodiment of the invention. As shown in FIG. 5, 

2 9 inverse transformer 560 includes frequency-to-time transformer 562 and overlap-and-add 

3 0 processor 564. Transformer 562 receives excised frequency domain vector XW* and transforms 
3 1 excised frequency domain vector XW' from the frequency domain to the time domain to produce 
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1 time domain vector "xw'." In one embodiment, transformer 562 uses an Inverse Fast Fourier 

2 Transform (IFFT) having a block size N according to the following Equation 5 : 
3 

4 (Equations) xw(n) = ^XW(k)e'^''^"''K, 

5 where 0<A:<(A^-1) 

6 XW(k) - frequency domain input; 

7 xw(k) - time domain output; and 

8 N= FFT block size. 
9 

10 Transformer 562 inputs time domain vector xw' to overlap-and-add processor 564. 

11 As shown in FIG. 5, overlap-and-add processor 564 receives time domain vector xw*. 

12 Overlap-and-add processor 564 adds a portion of a current time domain vector xw' to a previous 

13 time domain output vector to produce a current time domain output vector y*. In one 

14 embodiment, overlap-and-add processor 564 uses a "75% overlap-and-add process" to obtain 

15 current time domain output vector y*. The 75% overlap-and-add process adds the tail 75% of a 

1 6 previous time domain output vector y' to the head 75% of current time domain vector xw' to 

17 produce current time domain output vector y*. The current time domain output vector y' includes 

18 a current output block, which is the head 25% of current time domain output vector y\ Overlap- 

19 and-add processor 564 inputs current time domain output vector y* to output processor 170 of 

20 FIG. 1. 

21 Referring again to FIG. 1, output processor 170 receives time domain output vector y* 

22 from inverse transformer 160. Output processor 170 calculates valid time domain output vector 

23 y by multiplying the gain GW of windowed vector xw (FIG. 2) by the current output block. In 
2 4 one embodiment, gain GW equals 1 .68 for a Blackman window with 75% overlap-and-add 

2 5 processing. Output processor 1 70 outputs valid time domain output vector y. 
2 6 FIGURE 6 is a flowchart of an exemplary method that implements an embodiment of the 

2 7 invention. Certain details and features have been left out of flowchart 600 of FIG. 6 that are 

2 8 apparent to a person of ordinary skill in the art. For example, a step may consist of one or more 

2 9 sub-steps or may involve specialized equipment or materials, as known in the art. While STEPS 

30 610 through 670 shown in flowchart 600 are sufficient to describe one embodiment of the 
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1 present invention, other embodiments of the invention may utilize steps different from those 

2 shown in flowchart 600. 

3 As shown in FIG. 6, the exemplary method begins at STEP 610 in flowchart 600 whereat 

4 the method assembles time domain data blocks and transforms time domain data blocks into 

5 frequency domain vector XW. In one embodiment, the method uses an overlap-and-add process 

6 and a windowing technique to assemble time domain data blocks. In one embodiment, the 

7 method assembles time domain data blocks having a length N, wherein N is a power of 2. 

8 Exemplary overlap-and-add processes include 75% overlap-and-add processes. Exemplary 

9 windowing techniques include Blackman windowing according to Equation 1 . In one 

1 0 embodiment, the method uses a FFT to transform time domain data blocks into frequency 

11 domain vector XW according to Equation 2. After STEP 610, the method proceeds to STEP 

12 620. 

13 As shown in FIG. 6, at STEP 620, the method estimates power values of frequency bins 

14 of the frequency domain vector XW of STEP 610. In one embodiment, the method multipUes 

1 5 frequency domain vector XW by its complex conjugate to obtain power vector PXW comprising 

1 6 power values of frequency bins. After STEP 620, the method proceeds to STEP 630. 

17 At STEP 630, the method calculates a moving average of power vectors or average 

1 8 power vector PXW_avg from power vector PXW estimated at STEP 620 and selected previous 

19 power vectors PXW. According to the invention, STEP 630 scales each frequency bin of power 
2 0 vector PXW using power scale factor vector PSFV to obtain scaled power vector PXW*, which 
2 1 reduces suppression of high-ampHtude frequency bins of the current power vector when 

2 2 averaging power vectors. STEP 630 calculates average power vector PXW__avg from scaled 

2 3 power vector PXW*, selected previous power vectors PXW and the power scale factor vector 

2 4 PSFV. The method can discard the current power vector PXW from future average power vector 

2 5 PXW_avg calculations depending on transients. After STEP 630, the method proceeds to STEP 

26 640. 

2 7 FIGURE 7 is a flowchart of an exemplary method that implements an embodiment of 

2 8 STEP 630 of flowchart 600 of FIG. 6. Certain details and features have been left out of 

2 9 flowchart 700 of FIG. 7 that are apparent to a person of ordinary skill in the art. For example, a 

3 0 step may consist of one or more sub-steps or may involve specialized equipment or materials, as 
3 1 known in the art. While STEPS 732 through 740 shown in flowchart 700 are sufficient to 
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1 describe one embodiment of the present invention, other embodiments of the invention may 

2 utilize steps different from those shown in flowchart 700. 

3 As shown in FIG. 7, the method begins at STEP 732 whereat the method determines a 

4 power scale factor vector PSFV. In one embodiment of STEP 732, the method determines 

5 whether power values of frequency bins of the current power vector PXW are less than threshold 

6 T. An exemplary value for threshold T is shown in Equation 3. If a frequency bin power is less 

7 than threshold T, the method sets power scale factor PSF associated with the frequency bin to M, 

8 where M is the number of blocks used for averaging, else the method sets power scale factor PSF 

9 associated with the frequency bin to 1 . The method calculates a PSF for frequency bins of 

1 0 current power vector PXW to form power scale factor vector PSFV, which comprises a plurality 

11 of PSF. After STEP 732, the method proceeds to STEP 734 whereat the method calculates 

12 scaled power vector PXW by multiplying current power vector PXW by power scale factor 

13 vector PSFV. After STEP 734, the method proceeds to STEP 736. 

14 At STEP 736, the method calculates average power vector PXW_avg. In one 

1 5 embodiment, the method calculates average power vector PXW_avg using Equation 4. After 

1 6 STEP 736, the method can proceed to either STEP 640 of FIG. 6 or STEP 738. At STEP 738, 

17 the method determines current power vector PXW. In one embodiment, the method rescales 

1 8 scaled power vector PXW to current power vector PXW by dividing scaled power vector PXW 

1 9 by power scale factor vector PSFV. In one embodiment, the method obtains current power 

2 0 vector PXW from STEP 620 of flowchart 600. After STEP 738, the method proceeds to STEP 

2 1 740 whereat die method stores current power vector PXW along with previous power vectors. 

2 2 After STEP 740, the method proceeds to STEP 640 of FIG. 6. One of ordinary skill in the art 

2 3 shall recognize that the order of STEPS 736 and 738/740 can be switched without departmg from 

24 the scope and spirit of the present invention because current power vector PXW is calculated for 

2 5 use in calculating subsequent average power vectors PXW avg. For example, the method 

2 6 STEPS can have the following order: STEP 732, STEP 734, STEP 738, STEP 740 and STEP 

27 736. 

2 8 Referring again to FIG. 6, the method at STEP 640 excises selected frequency bins of 

2 9 frequency domain vector XW based on moving averages of bin values of average power vector 

3 0 PXW_avg. After STEP 640, the method proceeds to STEP 650. 

14 
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1 FIGURE 8 is a flowchart of an exemplary method that implements an embodiment of 

2 STEP 640 of flowchart 600 of FIG. 6. Certain details and features have been left out of 

3 flowchart 800 of FIG. 8 that are apparent to a person of ordinary skill in the art. For example, a 

4 step may consist of one or more sub-steps or may involve specialized equipment or materials, as 

5 known in the art. While STEPS 844 through 858 shown in flowchart 800 are sufficient to 

6 describe one embodiment of the present invention, other embodiments of the invention may 

7 utilize steps different from those shown in flowchart 800. 

8 As shown in FIG. 8, the method begins at STEP 844 whereat scaled median vector 

9 SMED is calculated. Scaled median vector SMED is calculated by multiplying median vector 

1 0 MED by median threshold MF T. Median vector MED comprises median values corresponding 

11 to frequency bins, where a median value for a specific frequency bin is determined by calculating 

12 the median of a set of frequency bins that includes the specific frequency bin. In one 

13 embodiment, the median value for a specific frequency bin is determined by calculating the 

1 4 median of a set of frequency bins comprising seven frequency bins tiiat are centered on the 

1 5 specific frequency bin. After STEP 844, the method proceeds to decision STEP 846 whereat the 

1 6 method determines whether an average power value of a frequency bin of average power vector 

1 7 PXW avg (i.e. , PXW_avg(n)) is greater than a scaled median value of a frequency bin of scaled 

18 median vector SMED (i.e., SMED(n)). If PXW_avg(n) is greater than SMED(n), then the 

19 method proceeds to decision STEP 850, else the method proceeds to STEP 848. 

20 At STEP 848, the method sets excision flag fExc(n) to 0. After STEP 848, the method 

2 1 proceeds to decision STEP 856. At decision STEP 850, the method determines whether previous 
2 2 excision flag Prev_fExc(n) is set to 1 . If previous excision flag Prev_fExc(n) is set to 1 , the 

2 3 method proceeds to STEP 852, else the method proceeds to STEP 854. At STEP 852, the 

2 4 method excises the frequency bin corresponding to PXW_avg(n) (i.e., set XW(n)) to 0. After 

2 5 STEP 852, the method proceeds to STEP 854. At STEP 854, the method sets excision flag 

2 6 fExc(n) to 1 . After STEP 854, the method proceeds to decision STEP 856 whereat the method 

27 determines whether excision flag vector fExcV calculation is complete. If not, the method 

2 8 returns to decision STEP 846. If so, the method proceeds to STEP 858 whereat the method 

2 9 outputs excision flag vector fExcV and excised frequency domain vector XW (i.e., frequency 

3 0 domain vector XW after excision). After STEP 858, the method proceeds to STEP 650 of FIG. 
31 6. 
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1 Referring again to FIGURE 6, the method at STEP 650 determines whether to include 

2 average power vector PXW in fixture average power vector PXW_avg determinations. In one 

3 embodiment, the method determines discard flag percent excised and eliminates transients. In 

4 one embodiment, the method determines excision percentage pExc from excision flag vector 

5 fExcV and determines whether current power vector PXW should be excised from future 

6 calculations of average power vector PXW_avg. In one embodiment, the method at STEP 650 

7 detects transient conditions based on a ratio between current excision percentage pExc and 

8 previous excision percentage pExc. After STEP 650, the method proceeds to STEP 660. 

9 FIGURE 9 is a flowchart of an exemplary method that implements an embodiment of 

10 STEP 650 of flowchart 600 of FIG. 6. Certain details and features have been left out of 

11 flowchart 900 of FIG. 9 that are apparent to a person of ordinary skill in the art. For example, a 

12 step may consist of one or more sub-steps or may involve specialized equipment or materials, as 

13 known in the art. While STEPS 952 through 966 shown in flowchart 900 are sufficient to 

1 4 describe one embodiment of the present invention, other embodiments of the invention may 

1 5 utilize steps different fi-om tfiose shown in flowchart 900. 

16 As shown in FIG. 9, the method begins at STEP 952 whereat the method calculates 

17 current excision percentage pExc. In one embodiment, the method calculates excision 

1 8 percentage pExc by dividing the summation of excision flags fExc of excision flag vector fExcV 

19 by length N. After STEP 952, the method proceeds to decision STEP 954 whereat the method 
2 0 determines whether current excision percentage pExc is greater than pExc_T, where pExc_T is 
21 the product of multiplying previous excision percentage pExc by a transient condition threshold. 
2 2 If pExc is greater than pExc_T, then the method proceeds to STEP 958, else the method proceeds 
23 to STEP 956. At STEP 956, the method sets exceed counter cExc to 0. After STEP 956, the 

2 4 method proceeds to STEP 966. At STEP 95 8, the method increments exceed counter cExc by 1 . 

2 5 After STEP 958, the method proceeds to decision STEP 960. 

26 At decision STEP 960, the method determines whether exceed counter cExc is less than a 

2 7 predetermined limit. In one embodiment, the predetermined hmit is three. If exceed counter 

2 8 cExc is less than the predetermined limit, then the method proceeds to STEP 962, else the 

2 9 method proceeds to STEP 964. At STEP 962, the method discards current power vector PXW 

3 0 fi-om fixture average power vector PXW_avg calculations. In one embodiment, the method 

3 1 discards current power vector PXW by setting discard flag F to 1 . After STEP 962, the method 
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1 proceeds to STEP 964 whereat the method sets current excision flag vector fExcV to previous 

2 excision flag vector Prev_fExcV. After STEP 964, the method proceeds to STEP 966 whereat 

3 the method proceeds to STEP 660 of FIG. 6. 

4 Referring again to FIG. 6, the method at STEP 660 performs a frequency domain to time 

5 domain transform on excised frequency domain vector XW. In one embodiment, the method 

6 uses an IFFT according to Equation 5. In one embodiment, the method uses an overlap-and-add 

7 process to assemble time domain data blocks. After STEP 660, the method proceeds to STEP 

8 670 whereat the method performs output processing and outputs data. In one embodiment, the 

9 method at STEP 670 calculates valid time domain output vector y by multiplying the gain GW of 

10 windowed vector xw by the current output block. In one embodiment, gain GW equals 1.68 for 

11 a Blackman window with 75% overlap-and-add processing. 

12 FIGURE 10 is a block diagram representing an exemplary ELF receiver application of a 

13 narrowband interference excision device according to one embodiment of the invention. As 

14 shown in FIG. 10, ELF application 10000 includes antenna 1010, pre-amphfier 1020 and 

15 receiver 1030. Antenna 1010 receives rf signals and is operatively coupled to pre-amphfier 

16 1020. Antenna 1010 inputs rf signals to pre-amplifier 1020. Pre-amplifier 1020 is operatively 

17 coupled to receiver 1030. Pre-amplifier 1020 inputs analog signals to receiver 1030. Receiver 

18 1030 includes analog-to-digital (A/D) converter 1032, narrowband interference excision device 

19 1034 and other receiver processing devices 1036. A/D converter 1032 is operatively coupled to 

20 pre-ampUfier 1020 and narrowband interference excision device 1034. A/D converter 1032 

21 receives analog signals from pre-ampUfier 1020, A/D converter 1032 converts analog signals to 

22 digital signals. A/D converter 1032 inputs digital signals to narrowband interference excision 

23 device 1034. Narrowband interference excision device 1034 is operatively coupled to other 

24 receiver processing devices 1036. Narrowband interference excision device 1034 removes 

2 5 narrowband interference from digital signals and inputs data to other receiver processing devices 

2 6 1036. In one embodiment, other receiver processing devices 1036 includes a clipper to suppress 

2 7 atmospheric noise. 

2 8 From the above description of the invention, it is manifest that various techniques can be 

2 9 used for implementing the concepts of the present invention without departing from its scope. 

3 0 Moreover, while the invention has been described with specific reference to certain 

3 1 embodiments, a person of ordinary skill in the art would recognize that changes can be made in 

17 



Express Mail No.: EE676641076US 



Navy Case No. 84436 



form and detail without departing from the spirit and the scope of the invention. The described 
embodiments are to be considered in all respects as illustrative and not restrictive. It should also 
be understood that the invention is not limited to the particular embodiments described herein, 
but is capable of many rearrangements, modifications, and substitutions without departing from 
the scope of the invention. 
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